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1. INTRODUCTION 


The primary objective of the ocean geoid Investigation was 
to determi-'e the fine-scale sea surface topography in the GEOS- 3 
calibration area using radar altimeter data. North-south and 
east-west components of the deflection of the v*rtlcal as well as 
values of the geoidal heights were estimated from an analysis of 
the dcta. 

The raw data files as provided by Wallops Flight Center 
(Ref. 1) contained all the necessary information to compute esti- 
mates of the geoidal heights from the GEOS-3 altimeter data. Only 
the altimeter data taken in the short-pulse mode were used for the 
investigation. Three major stages of processing were used in ob- 
taining the final results; the data flow diagram is shown in Fig. 
1. The first two stages used pass processors. Each satellite 
pass of altimeter data over the calibration area is individually 
processed by these programs. Two different filter processors were 
invv>stigated. The objective of the ti.iax processor was to combine 
all of the pass results to produce the geoidal heights and deflec- 
tions of the vertical. Two different methods were also investi- 
gated in the final stage of processing. 

A description of each processor is presented, as well as 
final results. 


Ref. 1. "GEOS-C Mission Plan," TK-63A0-01, Wallops Sta- 
tion, Wallops Island, VA. 


- 1 - 











WMtct d«l« 

K>« t*w |ri0 

— zn 


N t <l«f« 

•v«4w«t« fwnct«MH 

M fr*4 tctilwdi 
trcmrnfi 

^ r-* 


S*«M t >ut« 

wp*' 4Sd*(« 

tunclioni 

•f frxf lOAfiiudi 
ooennp 




A^«liv« g»o4l 

drNKfn>n«lHXi 
Junk«n« (•« «l I 

Mvthcxl |N«t« % 

«n<l •• 




Seyw»n d«ia 
•n<0 proc««<r^ 
btockt 

1 


r S 

f v«k*«t* 

Kinctkont «« 

tonpfudt 

croepofi 

^ J < ■ 


A / 

( 

tunctiOfH «t 
letitMl* 

I'OMinfi 


Fig. 1 Geoid Determination System Flow 





























2. DATA PREPROCESSOR 


The major objccCive of the data preprocessor was to Inter- 
face with the raw data as received from Wallops Flight Center and 
transform them Into a format more easily manipulated In the sub- 
sequent stages of processing. The raw data were used to compute 
estimates of the geoldal heights which were then aggregated and 
transformed Into measurements at equally spaced Increments of arc 
distance along the satellite subt'ark. A spacing of approximate- 
ly 6.4 km between measurements wai generated from the pass data, 
and data aggregation was achieved using a second-order polynomial 
smoother. The output of the preprocessor for each satellite pass 
consisted of equally spaced geoldal height measurementr. The lat- 
itude, longitude, and azimuth associated with each height mea- 
surement were also saved. The azimuth defined the satellite pass 
direction and was used by the final data processor. 
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3. FILTER PROCESSOR 


The object iv«i. 'f the second processor were to further 
snooth the output fron the preprocessor and to conpute the along- 
track deflection angles from the geoldal height measurements. 

Two different methods were Investigated: 

1. A Martin-Craham low-band pass filter (Kef. 2) and 

2. A filter based on the Wlener-Kolomogorof f theory 
(Ref. 3). 


3.1 MARTIN GRAHAM FILTER 

The Martin-Graham filter was the more flexible of the two 
methods considered. Since it gives a choice of cut-off frequen- 
cies, filter weights could be generated that would define the 
highest frequency passed by the filter. 


its 

ing 


The Martin-Graham filter is characterized by the fact that 
transfer function Is derivable from a transfer function hav- 

2 IT(W-W^) 

a roll-off of cos ■ where w Is the c»*.L-off frequency, 

w„ - w c 


the termination frequency. 


and w the Independent variable. 


TWo specific low-pass Martin-Graham filters were studied for use 
with the GEOS altimeter data. A smoother filter was used to gen- 
erate smoothed data on deflectlone of the vertical. In general, 
data having polynomial content, which Is not band-limited, would 
be effected In a strictly band-limited filter. To circumvent 
this undesirable effect, both filters were developed so as to 
pass without distribution linear functions of the height and 
slopes . 


Ref. 2. E. B. Anders, J. J. Johnson, A. D. Lasalne, P. W. 
Spikes, and J. T. Taylo, Digital Filters , NASA Report CR-136, 
December 1964. 

Ref. 3. A. Papoulls, Probability, Random Variables and 
Stochastic Processes , McGraw-Hill Book Co., New York, NY, 1965. 
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simulation atudlea were made to verify the effectiveneaa 
of the Martln*'Craham filters. The selection of the cut-off and 
termination frequencies Is of fundamental Importance for use In 
processing the altimeter data. These determinations ’’ere to be 
made by studying the spectral component of local fine-scale grav- 
ity surveys. However, prior to this It was decided to abandon 
this approach In favor rf the Wiener filter for reasons given In 
the following subsection. 


3.2 WIENER FILTER 

The method Investigated used a Ulener-Kolmogorof f nonrecur- 
slve filter based on closed covariance expressions derived by 
Tschemlng end Rapp (Ref. 4). Although not as flexible as the 
Martin-Graham filter. It was selected to process the data since 
(a) It was capable of filtering the data to the estimated noise 
level In the measurements (oi^ •• 24 cm) and (b) It preserved the 
physical significance associated with the covariance functions 
used In weight computation. 

Two sets of 35 weights were computed to process the data, 
one set for generating smoothed height data and the second for 
producing filtered deflection datri. The weights were generated 
so that the filtered result was based on data symmetrically 
placed relative to the filtered uolnt. Figures 2 and 3 are 
plots of the height and deflection weights, respectively, that 
were used. 

Using a typical pass of preprocessed altimeter data. Figs. 

4 through 6 Illustrate tiie type of results obtained with the Wiener 
filter. Figure 4 Is a plot of the gsoldal heights before filter- 
ing. Time Is used as the Independent variable rather than arc 
distance along the satellite subtrack. (One second Is equivalent 
to an arc distance of approximately 6.4 km.) The filtered heights 
are plotted In Fig. 5; both curves are shown on the same plot In 
Fig. 6; and the differences between the two heights are shown In 
Fig. 7. The filtered minus smoothed height residuals have a mean 
of -3.8 cm and sigma of 30 cm. Note the high-frequency structure 
In the residuals. The filtered along-track deflection angles are 
given In Fig. 8. 

The filtered height and deflection data computed from each 
satellite pass by the filter processor are archived for use by kl.c 
program used In the final stage. 


Ref. 4. C. C. Tschemlng and R. H. Rapp, "Closed Covarlan/*«* ^x- 
presslons for Gravity Anomalies, Geold Undulations, and Deflections 
of the Vertical Implied by Anomaly Degree Covariance Models," Ohio 
State University Report No. 208, May 1974. 
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Fig. 3 Wiener Filter Deflection Weights versus Tune. 
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Fig. 4 One Second Smoothed Geoidal Height Values (short pulse model). 
Latitude, longitude range from (30 . 299 ) to (41*’. 289 ). 
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Fig. 5 Filtered Geodial Height Values (Wiener filter, 35 weights, short pulse 
model). Latitude, longitude range from (30 , 299 ) to (41°, 289 ). 
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Fig. 8 Filtered Along Track Deflection Angles (Wiener filter, 35 weights, 
short pulse mode). Latitude, longitude range from (30 , 299 ) to 
(4r,289 ). 
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4. GEOID DETERMINATION PROCESSORS 


The objective of the final stage of processing was to com- 
bine all of the pass results to produce geoidal heights and both 
components of the deflection of the vertical. The data were pro- 
cessed using two different methods. Brief descriptions of each 
and the results obtained are given below. 


4.1 PROCESSOR I 

It was Intended that the first processor be made as simple 
as possible, yet provide reasonable value checks for ‘ \e subse- 
quent processors. Thus the following approach was implemented. 


4.1.1 Data Selection 

Filtered passes outside the calibration grid were eliminat- 
ed. Within each puss only data confined to the calibration grid 
were kept. Following the selection process, the pass boundaries 
were Ignored. Processor I had two additional functions at this 
stage: 

1. To separate south- and north-going data and 

2. To evaluate the functions at latitude or longitude 
grid intersections. 

The functions evaluated were geoidal height, deflection angle, and 
either latitude or longitude. The "evaluation" was a simple aver- 
aging of all data in the vicinity (+1*) of each latitude or longi- 
tude grid intersection. Noisy data were eliminated using 3-slgma 
testing. 


4.1.2 Sorting 

A standard sorting program was used that took the data and 
sorted them into latitude, longitude or longitude, latitude order. 
With Processor I, the data have already been evaluated at either 
latitude or longitude grid intersections. However, for Proces- 
sor II, all the data must be sorted. 
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4.1.3 Function Evaluation 

The next nodule In Procesaor I operated on Che sorted data 
Co evaluate the functions at either latitude or longitude grid 
Interaectlons . The evaluation was a simple averaging of all data 
In the vicinity (fl*) of each latitude or longitude grid Intersec- 
tion. (If data are at latitude grid Intersection, then the aver- 
aging la for longitude grid crossing and vice versa). The output 
of this program Is a table (file) of geoldal heights and deflec- 
tions at grid Intersections. 


4.1.4 Generation of Geoidal Map' 

The final step In Processor I Is to combine the four Cables 
(files) and to generate a map (table) of geoldal heights and east- 
west and north-south deflections of the vertical. The north-south 
pairs are averaged and then combined to generate the east-west and 
north-south deflection angles. The geoldal height Cable Is produced 
by averaging the north and south geoldal heights. The results ob- 
tained with this processor are shown In Figs. 9 through 11. 


4.2 PROCESSOR II 

The second method Incorporates a wc'..htlng function approach 
for the modeling of Irregular surfaces developed by Junklns et al. 
(Refs. S and 6). This method has many f>>itures desirable In pro- 
cessing large quantities of data. Tt al lo«? the user to partition 
Che data Into a sequence of overlapping ai.usets, each of which Is 
processed separately. A separate model Is used In fitting each 
subset of data Chat Is applicable only for that subset, thus re- 
quiring overall a less sophisticated model. By Che use of weight- 
ing functions, nth-order continuity across boundaries of Che sub- 
sets can be enforced. A complete explanation of the method Is 
given In Refs. 5 and 6. 


Ref. 5. J. L. Junklns. G. W. Miller, and J. R. JancalCls, 

"a Weighting Function Approach to Modeling of Irregular Surfaces,' 
J. Geophys. Res. , Vol. 78, No. 11, 10 Apr 1973. 

Ref. 6. J. R. .Tancaltls and J. L. Junklns, "Modeling In n 
Dimensions Using a Weighting Function Approach," J. Geophys. Res. , 
Vol. 79, No. 23, 10 Aug 1974. 
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4.2.1 Segmentation of Data Into Blocks 

Processor II requires as Input a sorted « ata fil#: (latitude, 
longitude order) where all data are restricted to the calibration 
grid. An output file Is generated where data have been segmented 
Into grid areas with each group of data (grid block) preceded by 
a special header record that Identifies the block (beginning and 
extent of latitude and longitude) and gives a count of the number 
of data points (records) In the block. The output file (segmented 
data file) Is the Input for Processor II. 


4.2.2 Function Evaluation A 

In each region, the filtered height measurements are repre- 
sented as nth-order polynomials 

"i-Z ecv 

t»0 m*0 


where 


■ <p - <Pj , 

■ X - , 

9 ( 1 )X( 1 )- reference origin associated with the 1th region, 
and 

fit parameters arsoclated with the 1th region, 
tn 

Ihe parameters are determined from a least-squares fit to the 
filtered measurements. 

Components of the deflection of the vertical are defined ca 
1 dH 

^ - — (n-s component), (2a) 

r cfp 

n ■ - ^ (e-w component), (2b) 

r cos 9 a 


In tft and X : 

AXj"’ , (1) 
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Differentiating Eq. 1 yields 


L-1 L-l-1 

-tL E * 

t*0 m«0 


I-*-! ,m 


(t+15 AX 


(3a) 


L-1 L-l-1 


Y' y* A, (m^l) AX® . (3b) 

r coa tfi t,ort-l ^ 


l*C p«0 


Thus the fitted values are also used to define the ( and n 

components . 

The results of this determination are presented In Figs. 
12 through 14. 


4.2.3 Function Evaluation B 


This final ivcthod also uses the weighting function algo- 
rithm described in Refs. 5 and 6. However, here the A. para- 
meters are determined from a least-squares fit to the filtered 
along-track deflection angle measurements. In each region, the 
deflection angle is modeled as 


A. ■ -(C cos A n sin A ) , (4) 

X z z 

where A Is the azimuth angle associated with the measurement 
point. ^The filtered along-track deflection angles are computed 
such that 



(5) 


where s defines distance along the satellite subtrack (thus, the 
need for the minus sign In Eq. 4). The functional forms for C and 
n are given by Eq. 3. 
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In chl« least-squares fit, the A_ parameter Is not esti- 
mated. Thus, the height as defined In Eq. 1 can only be computed 

minus the constant^ (l.e., ■ H . 

A value for Is then estimated by computing the mean 

( 1 ) 

value of H - AH for the region being processed, 

in 

The results of this method are preaented In Figs. IS through 
17. 
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5. CONCLUSIONS 


The results obtained with the three different methods are 
very similar. The rms height difference among the three methods 
is only 0.5 m. The results also compare favorably with gravi- 
metrically determined heights in the calibration area as well as 
with those obtained by Hadgigeorge and Trotter (Ref. 7). The rma 
of the deflection differences between the last two methods were 2 
ai.'d 1.2 arc sec in C and n« respectively, versus 3 and 2.5 arc 
sec for the differences with method one. 

In the final two methods, second-degree polynomials in Sf 
and X were used to model the height in each region. The rms of 
the fits to the filtered height and the along-track deflection data 
were 3.1 m and 3.2 arc sec, respectively, with both methods. The 
use of a higher degree polynomial in the fit did not significantly 
improve these results. 


Ref. 7. G. Hadgigeorge and J. E. Trotter, "Short Arc Reduc- 
tlcns of GEOS-3 Altimetric Data," Geophys. Res. Lett. , Vol. 4, 

No. 6, Jun 1977. 
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